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DAISY - AN EFFICIENT PROGRAM FOR PRODUCTION OF NEUTRON 
PSEUDO RESONANCE PARAMETER CHAINS IN THE 
UNRESOLVED ENERGY REGION 
by Thor T. Semler 
Lewis Research Center 

SUMMARY 

The code DAISY provides a fast technique for generating pseudo resonance param- 
eters in the unresolved resonance region. The code makes use of the Central Limit 
Theorem in order to select chains of resonance parameters whose average values are 
inside a predetermined limit of accuracy. The Kolmogorov- Smirnov test for goodness 

of fit is used on chains of 5 through 100 resonances to determine the chain that best fits 

o 

the Porter-Thomas and Wigner distribution functions. The x test for goodness of fit 
is used on chains of more than 100 resonances to determine that chain which best fits 
these distribution functions. 

The code DAISY is found to be at least an order of magnitude faster than PSEUDO, 
a comparable code. 

As an example, the code DAISY is used to compute s- and p-wave resonance param- 
238 

eters for U near 20 keV. From these parameters continuous cross sections are 

9 QO 

computed and used to calculate the Doppler effect in a 0. 0155-centimeter-thick U 

238 

foil and in a 2-centimeter-diameter U sphere, at 20 keV. 

It is found that using the set of pseudo resonance parameters chosen by DAISY, one 
may preserve the measured value of the total cross section, and the Doppler effects so 
computed, are compatible with other theoretical computations of the Doppler effect. It 
is also found that, while the average s-wave capture cross section at 20 keV is about 
one-third that of the p-wave, the s-wave contribution to the Doppler effect is about 
60 percent. The effect of neglecting resonance series overlap is to overestimate the 
Doppler effect. 



INTRODUCTION 


Continuous neutron cross sections in the unresolved energy region are needed to 
calculate fast neutron Doppler effects, to analyze fast neutron transmission experiments, 
and to estimate stellar reaction rates. The unresolved resonance region generally spans 
an energy region from about 1 keV to about 1 MeV. In figure 1 the neutron capture cross 
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Figure 1. - Portion of experimental capture data from first U 238 sample in neutron beam from Petrel event. Small arrows indicate loca- 
tions of weak resonances which are assigned 1=1. Data from more sensitive detectors and thick sample were used in confirmation 
and measurement of these resonances. 


238 

section of uranium-238 (U ) is shown in a part of the resolved resonance region 

238 

(ref. 1). In figure 2 the neutron capture cross section of U is shown for the unre- 
solved resonance region (ref. 2). Figure 2 depicts the results of 16 different experi- 
ments that differ by as much as 40 percent. The results of figure 2 are averages over 
several resonances. 

Various analytic approaches have been adopted to calculate average cross sections 
in the unresolved region. In an attempt to evaluate fission product reactivity effects, in 
intermediate reactors, a method using strength function data was developed (refs. 3 
to 5). However, this statistical averaging method does not provide the detailed energy 
dependent cross sections necessary for the calculation of a fast reactor Doppler effect. 

If it is assumed that the statistical distributions of resonance parameters found at 
low neutron energies are obeyed at the higher unresolved energies, it is possible to 
sample from these distributions and generate "pseudo” or "mock" resonance param- 
eters. These pseudo resonance parameters have been used to compute energy dependent 
cross sections with correct average properties (refs. 6 to 10). Brissenden and 

238 

Durston have used sampling techniques to generate sets of mock cross sections for U 
(ref. 6). 

238 

In a series of experiments designed to measure the Doppler effect in U , Perkin 
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Energy of neutron, E N , MeV 

Figure 2. - Portion of unresolved resonance region for captures. 


and Fieldhouse used a modification of the method of Brissenden and Durston to compute 
the Doppler effect theoretically (refs. 11 and 12). However, the computer time neces- 
sary to generate these mock cross sections is prohibitively large, ~2. 5 hours per iso- 
tope per temperature on an IBM- STRETCH computer. 

The code DAISY described herein is one to two orders of magnitude faster in the 
generation of resonance parameters than the code PSEUDO (ref. 9) which has a timing 
of about 139 resonances per second on the UNIVAC 1108. DAISY also allows substantial 
flexibility, in that, short resonance chains, as few as five resonances in a chain, can 
be generated. 

In order to illustrate the versatility of the code DAISY, it is used with the code 
GAROL (ref. 13) to analyze the effect of s- and p-wave resonance series overlap. The 
method is applied to compute the Doppler effect for a 0. 0155-centimeter slab of 
and a 2. 0-centimeter sphere of at 20 keV (refs. 12 and 14). 


STOCHASTIC GENERATION OF RESONANCE PARAMETERS 


In order to generate neutron resonance cross sections, one needs the resonance pa- 
rameters associated with the Breit-Wigner single-level formula. These are the energy 
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of the resonance E Q , the scattering width of the resonance r n , and the capture width 

of the resonance r . 

y 

These measured parameters allow one to compute accurate values of the continuous 
neutron cross sections from resonance to resonance (ref. 13). 

It has been asserted on the basis of theoretical considerations that these resonance 
parameters have definite analytic probability distributions (refs. 15 to 17). The level 
spacing distribution of Wigner and level width distribution of Porter and Thomas have 
been checked experimentally in the resolved resonance region for many isotopes, and 
they have been found to adequately describe the distributions of parameters measured 
experimentally (refs. 18 to 20). 

One may assume that the statistical distributions found to obtain at low neutron 
energies are also valid for higher neutron energies. Thus one may employ these dis- 
tributions of resonance parameters to generate a set of pseudo resonance parameters 
in the unresolved resonance region. The code DAISY described below allows one to gen- 
erate chains of pseudo resonance parameters, in the unresolved resonance region, 
which preserve the correct average values, as well as, an accurate representation of 
the correct distribution functions. 


Wigner Level Spacing Distribution 

It is found that the level spacing between two neutron resonances, of like J and Z , 
follows a probability distribution of the form 

p(x) = - x exp (— x^ x > 0 
2 V4 J 

p(x) = 0 x < 0 (1) 

where J is the spin quantum number of the compound nucleus formed by the target and 
the neutron, Z is the angular momentum of the neutron with respect to the target nu- 
cleus, and x is a normalized dimensionless level spacing or 

x = J- (2) 

D J 

where S is the nearest neighbor level spacing in eV and Dj is the J state average 
level spacing in eV. Equation (1) is then the Wigner level spacing distribution (refs. 

16 and 17). 
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Level Width Distribution of Porter-Thomas 


The variation of the reduced scattering widths r° is found to follow a x^ distri- 
bution function (ref. 18) with one degree of freedom for s-wave ( l - 0) resonances while 
two degrees of freedom have been suggested for p-wave (l = 1) resonances (ref. 15); 
that is, 


■pO = Ip 

r n E 0 v l r n 


(3) 


where E Q is the resonance energy in eV, ^ is the penetrability factor, and T n is the 
scattering width associated with the Breit-Wigner formula. The distribution of 


y 



(4) 


is found to be of the form 


p(y, v) 


2^/ 2 r p! jj 1 y( t '/2)- 1 e -y/2 



p(y, v) = o y < 0 


(5) 


where the number of degrees of freedom v is equal to 1 for s-wave resonances and 
equal to 2 for p-wave resonances. This is the distribution of Porter and Thomas for 
scattering widths. 

The capture width is found to vary little from resonance to resonance. This, 

in fact, is consistent with the theory of Porter and Thomas which indicates r should 

2 ^ 
be distributed as x with a large number of degrees of freedom (ref. 21). Since the 

o 

mean square deviation about the mean for a normalized x distribution is two over the 
number of degrees of freedom, r 1 should remain relatively constant from resonance to 
resonance. However, some recent experiments tend to cast doubt on the assumption of 
a constant (refs. 1 and 22). 

If these distributions are assumed to be valid for higher neutron energies, one may 

sample from the Wigner level spacing distribution to obtain a series of E n f s and then 

2 ^ 
sample from a x distribution (Porter-Thomas) for the corresponding r n T s. This 

technique may be used to generate continuous cross sections in the energy region be- 
tween the highest measured resonance and the energy region where the Doppler width 
becomes larger than the level spacing. 
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Sampling From the Normalized Wigner Distribution 


When the Wigner level spacing distribution is written as equation (1), the Dj mul- 
tiplier has been normalized out. The expectation value of x, or the mean of x, is 1. 
It is adequate, for the present, to discuss the sampling from this normalized probabil- 
ity function. In order to sample X, a level spacing, from the normalized Wigner dis- 
tribution (eq. (1)), one forms the functional mapping (eq. (6)) 


co 


p(x)dx 


( 6 ) 


where is a random number uniformly distributed on the interval (0, 1) 


R "i/ Xxexp (i x2 ) dx 


(7) 


Solving the previous equation for X gives 


X = A = J-ln(l - R n ) 

V7 T 1 


( 8 ) 


or since 1 - is distributed as R n 


X = J- ^-ln(Rj (9) 

V7T 

One may now easily gain the correct J state level spacing by multiplying X by Dj, 
the J state average level spacing. 

Sampling From the x 2 or Porter-Thomas Distribution 

The Porter-Thomas or neutron width distributions can likewise be written in a 
normalized form 


p(y, 1) = jV ^ 2 Jj y " 1 /2 e -y/2 (s-wave) 


( 10 ) 
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(p-wave) 


(ID 


p(y, 2) = [2r(l)] _1 2 _1 e“ y/2 

where y is given by equation (4). It can be shown (see appendix A) that the problem of 

o 

sampling a member from a x distribution with v degrees of freedom can be reduced 
to the problem of sampling u random values selected from a normal distribution of 
mean 0 and variance 1. That is, if the sequence 


Y-|, Y 2, Yg, . . . , Y n , . . . 


is composed of values selected from a normal distribution of unit variance and mean 0, 
the sequence of transformed random variables 




has a 


distribution with one degree of freedom and the sequence 


( y i + y !)> + • • •• ( y L-i + y L) • • ■ • 

2 

has a x distribution with two degrees of freedom. However, each element of the sec- 
ond sequence must be divided by two in order to preserve normalized property. There 
are available many techniques for sampling from the normal distribution. A fast algo- 
rithm with a correct treatment of the wings of the normal distribution has been chosen 
(ref. 23). It is described in appendix B. 

Once a series of Eq's or resonance energies have been generated, one can sample 
o 

from the appropriate x distribution; the value of T n associated with an Eq is then 

r n = E 0 /: Vn Y n (s-wave) (12) 


or 

r n - E 0 /2l, l Y n ( Y ln- 1 + Y 2n) rl <P- wave > 
where vq and are the penetrability factors, 

uq - 1 (s-wave) 


(13) 


(14) 
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v-i = (p-wave) (15) 

1 + k 2 R 2 

R is the nuclear radius and k is 2jt times the wave number of the incident neutron 
(ref. 24). 

Although it is possible to produce such resonance chains at will, a question arises 
as to how one might compare one such chain to another. Of course, given the same set 
of initial average resonance parameters, one will generate a somewhat different chain 
each time a chain is generated, depending upon the sampling method. For example, 
given a different random number generator for sampling purposes, one might use the 
same technique twice and obtain different results. Stated more accurately, given n 
such generated chains of m resonances, which of the n chains should be chosen as a 
representation of the actual resonance structure in this energy region? This difficulty 
has sometimes been rejected, by asserting, that since one was sampling correctly from 
the distributions, one chain of resonance parameters, so generated, was as good as any 
other such chain (ref. 6, pp. 75-76). Alternatively, an integral approach has been 
taken (ref. 9). In one such integral approach, several resonance chains are calculated, 
the (dilute resonance) integral of each chain is computed, and a chain with a resonance 
integral near the mean of these values is chosen as the best chain. 

A different approach has been adopted herein, which gives not only the correct reso- 
nance integral value but also provides a much faster running time - comparable with 
single chain times. The present method for the choice of the best chain depends upon 
the use of the Central Limit Theorem, and on two statistical techniques for goodness of 
fit that are discussed separately. One technique has been devised for very short chains 
of resonances, 5 to 100 spacings. The other technique is applied to longer chains of 
101 to 5000 spacings. 

Central Limit Theorem 

The expectation value of a sum is the sum of the expectation values, that is, if 
gj, gg, . . . , g n are any random variables which have expectation values, then 

E(g x + g 2 + • • • + g n ) = E(gj) + E(g 2 ) + . . . + E(g n ) (16) 

If each of g. are independently sampled from the same distribution, then the following 
relation results: 

E(g x + g 2 + • • • + g n ) = nE(gj) = nE(g 2 ) = . . . = nE(g n ) (17) 
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Also if gp gg, . . . , g n are independent random variables, the variance of the sum is 
equal to the sum of the variances. That is, 


gf+g 2 + 


+g 


= cr“ + cr 


n 


H §2 


+ cr 


3 n 


(18) 


Again, if each of the g^ are independently sampled from the same distribution, 


gl+g 2 + 


+g 


= no 


n 


Si 




(19) 


The Central Limit Theorem (refs. 27 and 30) states that, if gj, gg, . . . , g n are 

random variables independently chosen from the same distribution each with expectation 

2 

p and variance a , then the probability distribution of the sum 


S n = gi + g 2 + 


• + g 


n 


is asymptotically normal as n — °° or 


lim P 
n~°° 





-t* 72 


dt s N(w) 


( 20 ) 


( 21 ) 


It is found empirically that this limiting distribution, the standard normal distribution is 
assumed very rapidly for X^'s chosen from the Wigner distribution. In fact, for n as 
small as 5, the distribution of the sums of X^’s chosen from the Wigner distribution 
provide a good approximation to the standardized normal distribution. Thus one can 
compute both the probability, that the sum of n samplings from the Wigner distribution 
will be below a value S_ as well as the probability that the sum of n samplings will be 
above some value S . This calculation allows one to select a probability level, say 
85 percent, and to reject a chain immediately if it does not approximate the average 
level spacing better than 85 percent of such chains. At this point, it is useful to work 
through an example of such limits. The initial conditions are that 

E initial = 4500 eV 


E final 


= 5000 eV 


Dj = 5 eV 
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Since 


E final 



initial 


100 


the distribution function of the sum of 100 samplings and its variance from the normal- 
ized Wigner distribution is evaluated. It is well known that the variance is the expecta- 
tion of the square of the difference between the mean and the random variable or 


CT 2 = E[(X - ju) 2 ] 

(23) 

a 2 = E(X 2 ) - , u 2 

(24) 


Thus for the normalized Wigner distribution 


< r 2 = f x 2 p(x)dx - l 2 

•'a 


(25) 


•'O 


x 2 - x exp^— x 2 ^dx - 1 


(26) 


a 2 1 = 0.273239 . . . 


(27) 


<r = 0. 5227232 


(28) 


One can now generate the standard normal distribution, which very closely approximates 
the distribution of the sum of 100 random variables S^qq sampled from the Wigner dis- 
tribution after the form of equation (21); that is, 


S 100 " 100 x 1 \ 

< w) 

0. 5227 v/100 / 


i N(w) 


P 


S 100 " 100 \ 
< w 

5.227 / 


N(w) 
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One can now answer such questions as what is the probability that the sum S^qq is be- 
tween 99 and 101; that is, 


g /99 - 100 ^ _ 0 19 \ _ N (_ 0 19 ^ _ 0 4247 
\ 5. 227 J 


P( 101 r as 0. 19) = N(0. 19) = 0. 5753 
\ 5. 227 j 

P(99 < S 100 < 101) = 0. 5753 - 0. 4247 = 0. 1506 

If, for example, a probability level of 85 percent had been chosen for a chain of 100 
levels, any chain adding up to less than 99 or more than 101 could be ignored. More de- 
tailed calculations could then be used to choose the best chain from the 15 percent ac- 
cepted by this technique. Also, since the normal distribution is a symmetric function 
about the expected value of the sum, no bias has been introduced. 

A similar technique is used for the scattering widths. Since the level widths are 
2 2 

distributed as x > and as x does not approximate a normal distribution until many 

degrees of freedom are introduced, the sampling limits are chosen from the appropriate 
2 

X distribution. Again the expectation of the sum is equal to the sum of the expectations. 
The sampling limits are adjusted such that no bias is introduced; (i. e. , equal probability 
about the expected value). For an example similar to the aforementioned, the expected 
value of the sum of 100 values of is 100 and the 85 percent limits - 15 percent 

acceptance - are 97. 375 and 102. 715. The unequal spacing about 100 of the limits is 
such that no bias is introduced into the accepted values (ref. 26). 

The aforementioned technique allows one to generate many resonance chains, all 
drawn from the correct sampling distributions, and to immediately reject those chains 
which are found to vary from the average by more than the limits set. No extensive 
calculations are wasted on these rejected chains. The values of both D T and r° are 
preserved to within set limits. In the next section two techniques will be outlined which 
allow one to choose from the accepted chains, that chain, which best fits the distribution 
functions for D and r . 

FINAL SELECTION OF BEST RESONANCE CHAIN 

In the preceding section it was shown how one can reject resonance chains prior to 
any involved or detailed considerations on strictly statistical grounds. Two methods 
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which allow one to choose, from chains found acceptable by the above test, the "best" 
chain, in the sense of goodness of fit, are now described. The method found applicable 
to short chains (where n, the number of resonances, is less than or equal to 100 but 
greater than or equal to 5) was the Kolmogorov-Smirnov statistic for goodness of fit. 


Kolmogorov-Smirnov Test 

The cumulative distribution function F(x) (eq. (29)) equals the probability that the 
random variable has some value less than or equal to x; that is, 


F(x) = f p(x)dx 

•<_ on 


(29) 


If one also forms the discrete cumulative distribution function F n (x) for n samples 
from the distribution function, one can compute the Kolmogorov-Smirnov statistic A 
(ref. 27) as 


A = least upper bound of | F n (x) - F(x) | - 00 < x < 00 (30) 

The meaning of this statistic is illustrated in figure 3. Once one has computed A for 
each acceptable chain, that chain with A a minimum is chosen as the "best" chain. It 
can be seen that the minimum Kolmogorov-Smirnov statistic corresponds to the best 



Figure 3. - Illustration of Kolmogorov-Smirnov test for goodness of fit. 
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fitting of the continuum distribution function by a discrete function. Thus by the intro- 
duction of the Central Limit Theorem the measured mean values of the level spacing and 
the scattering width are conserved, and by use of the minimum Kolmogorov- Smirnov 
statistic a very good fit to the distribution function is obtained. Hence, all the informa- 
tion available to the construction of the "best" chain has been applied. This method is 
quite fast for chains of n resonances where 5 < n < 100. 

For chains containing more than 100 resonances, a different technique has been em- 
ployed in order to facilitate the choice of the "best" chain. The method which was 

o 

adopted for the longer chains is the x test for goodness of fit. 

X 2 Test for Goodness of Fit 

Due to the necessity of sorting individual values to compute F n (x) for the aforemen- 
tioned Kolmogorov- Smirnov test, it is not particularly well suited to the production of 
long resonance chains. Again, however, we should like to maintain the correct distri- 
bution function. The technique best suited to such a requirement for large samples is 
2 

the x test for goodness of fit. For this particular application, it has been imple- 
mented in the following manner. The probability distribution has been divided into 20 
equally probable portions. Twenty were chosen in order to maintain an expected value 
in each portion, or bin, of more than five (ref. 28) Once a chain has been accepted by 
the Central Limit Theorem, the random variables used to generate it are placed in the 
appropriate bin corresponding to the 20 equally probable ones. Finally, when all the 
random variables have been tallied, the expected value for each bin is subtracted from 
the number actually found. The twenty values so generated are squared and summed. 
This sum is then divided by twenty and multiplied by the number of resonances. The 

resultant number is called the x^ °f the fit and is a measure of the goodness of fit. If 

2 . 2 
X is large, the fit to the distribution is poor; likewise, if x is small, the fit is good. 

Hence, the chain with minimum \ is chosen as the best fit to the distribution function. 

This x 2 is a number and is not to be confused with the distribution function of x 2 which 

was discussed earlier in this report. The connection of this number and the distribution 

function is discussed in most books on statistical methods (refs. 29, 30, 25, and 28). 

Hence, the DAISY program by the use of rather simple and fast techniques has chosen 

the resonance chain which has an average level spacing and average reduced scattering 

width which are inside predetermined limits. DAISY has also chosen that chain whose 

elements best fit the two distribution functions associated with level width and level 

spacing. 

Figure 4 shows the mechanics of the code DAISY in block diagram form. 
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Figure 4. - DAISY flowchart. 


Timing 

In summary, the code DAISY provides a technique for the simultaneous production 
of correct average quantities and accurate statistical distributions in a short period of 
computing time. 

As an example, 1000 chains of 58 p-wave resonances with a rejection probability 
of 95 percent were generated by the code DAISY and the best chain chosen in a total 
computer time of 0. 71 minute on an IBM 7094-11. This gives a speed of 1362 resonances 
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per second on the IBM 7094-11. This may be compared with the timing of the PSEUDO 
code (ref. 9) which generates 1248 resonances in 9 seconds on the UNIVAC 1108 or about 
138. 7 resonances per second. Since the UNIVAC 1108 has a cycle time of about 0. 75 
microsecond and the 7094-11 has a cycle time of about 1. 4 microsecond, one may make 
a rather crude comparison of the two machines by the ratio of cycle times. If this is 
done for the PSEUDO code, one obtains a value of about 74. 2 resonances per second on 
the IBM 7094-II. On the basis of this admittedly crude comparison, the DAISY technique 
is at least, since p-wave calculations are longer, 18.4 times faster than the technique 
of PSEUDO. 

In comparing problem size limitations the maximum number of chains for DAISY is 
9999, and the maximum number of ladders for PSEUDO is 500. The maximum number 
of resonances per chain in DAISY is 5000, and the maximum number of resonances per 
ladder for PSEUDO is 4000. 


CALCULATION OF U 238 RESONANCE PARAMETERS NEAR 20 keV 

The use of the DAISY code for the generation of resonance parameters is illustrated 
238 238 

for U . U has both zero spin and relatively well measured s- and p-wave strength 

functions (refs. 1, 20, 31 to 34). The zero spin property allows one to sample from a 

single Wigner level spacing distribution for s-wave resonances. The neutron energy 

region about 20 keV was chosen since the contribution of the p-wave component to the 
238 239 

U (n, y) U cross section is near its maximum and should allow a convenient check 
on both the s- and p-wave components of the cross section averages. 

To provide a set of resonance parameters suitable for further calculations, the fol- 
lowing procedure is used: 

(1) Using the strength functions and the average level spacing, r° (s-wave) and r° 
(p-wave) are determined. 

(2) Within the energy limits, Ej n ^ a j to Eji na p the sets of resonance parameters 
r n (s-wave) and Eq (s-wave) are generated by DAISY for the best chain (I = 0) or 
chains (I # 0). 

(3) The measured value of and the values of r° (p-wave) and Dj (p-wave) 

are used to generate the sets of resonance parameters T n (p-wave) and Eq (p-wave) in 
this energy region for the best chain (unseparated spin statistics) or chains (separated 
p-wave spin statistics). 

(4) Using the resonance parameters just generated and Op 0t (measured), average 
partial cross sections are computed for the range E initial to E final - 
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(5) The sum of the average partial cross sections 


a tot( computed ) = °c + CT c + ^s 


+ a g + Op ot (measured) 


is compared with a^ o ^.(measured) in this energy region. 

(6) If necessary, the measured value of a pot is adjusted such that u tot (measured) 
equals ^(computed). 

(7) The final set of resonance parameters (p-wave parameters must be changed due 
to eq. (15)) will match the average total cross section and will have the specified aver- 
age parameter values within the set limits. 


RESULTS 

Average Cross Section Calculations 

O QQ 

As a test of the technique used in DAISY, the average partial cross section of U 
are calculated by the code GAROL (ref. 13). The GAROL code assumes a lattice of 
absorber lumps embedded in a moderating medium. It solves two coupled integral 
equations for the fluxes in each region. It also allows one to compute Doppler broadened 
cross sections directly from the Breit-Wigner resonance parameters provided by DAISY. 
A hydrogen- like (A = 1.0) pure scatter er has been chosen for the moderating material. 

A narrow energy region, 19. 5 to 20. 5 keV, has been chosen. 

DAISY is used to generate s- and p-wave resonance parameters from 19. 5 to 20. 5 
keV. The information used by DAISY is shown in table I. The p-wave parameters have 
been generated in this energy region using a measured value of cr pot equal to 10. 8 barns 
(ref. 35). The p-wave parameters have been chosen as merged, no J dependence, as 
there is no measured data on separated spin statistics for the p-wave resonances of 


TABLE I. - PARAMETERS USED IN DAISY FOR U 238 


Parameters 

s-Wave 

p-Wave 

Angular momentum of neutron with respect to target nucleus, l 

0 

1 

Capture width, r , meV 

23. 8 

23.8 

Average reduced scattering width, meV 

1. 83 

3.7 

Average level spacing, D, eV 

20. 8 

7.0 

Potential scattering cross section, °p ot , b 

10.8 

— 
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Figure?. -S-wave capture cross section. 
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Figure 6. - S-wave total and potential cross section. 


The s-wave capture component and the s-wave total cross section showing the scat- 
tering interference with the potential cross section about 20 keV are shown in figures 
5 and 6. A portion of the p-wave capture cross section and the p-wave total cross sec- 
tion are shown in figures 7 and 8. The total cross section and the total capture cross 
section are shown in figures 9 and 10. These figures (figs. 5 to 10) show cross sections 
which have been Doppler broadened to 300 K and have been computed at 0. 5-eV inter- 
vals. 

Table II shows the average values of the cross section components of the average 
total cross section at 20 keV for an initially selected value of 10. 8 barns for the poten- 
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Neutron energy, eV 

Figured - Portion of p-wave capture cross section. 



Neutron energy, eV 

Figure 8. - Portion of p-wave total cross section. 
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Figure 9. - Total cross section. 



Figure 10. - Total capture cross section. 


tial cross section. In column 3 of table II the results of an analytic computation of these 
averages is shown (ref. 36). This analytic method, that of Lane and Lynn, allows one 
to accurately compute average cross sections, both s- and p-wave. The results of the 
DAISY-GAROL computation are in quite good agreement with these analytic results. 

This agreement indicates the sampling techniques used in DAISY have been implemented 
correctly. However, according to steps (5) and (6) in the previous section, it is neces- 
sary to choose a value of a so as to maintain the measured value of The cross 
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TABLE n. - COMPONENTS OF AVERAGE TOTAL CROSS SECTION FOR U 238 AT 


20 keV BY THE METHOD OF DAISY AND THAT OF LANE AND LYNN (REF. 36) 


Component 

DAISY 

Lane and Lynn 
(ref. 36) 

DAISY 

(FINAL) 

Potential scattering cross section, 0 ^^, b 

10. 800 

10. 800 

9.647 

Average s-wave scattering cross section, or , b 

s 

2.378 

2.381 

2.380 

Average s-wave capture cross section, a , b 

s 

Average p-wave scattering cross section, a , b 

. 165 

. 163 

. 165 

1.226 

1. 223 

i 

1. 074 

— P 

Average p-wave capture cross section, a , b 

. 536 

. 527 

. 514 

— P 

Average total capture cross section, a c , b 

.701 

.690 

.679 

Average total cross section, v^ot’ b 

15. 105 

15. 094 

13. 780 


sections so generated are shown in the last column labeled DAISY (FINAL). 

Table III shows three values of the measured average total cross section at 20 keV 
and three values of the potential cross section inferred as necessary to maintain these 
three total cross sections (see step (6) in previous section). As it would be extremely 
difficult to assign any one of these measurements a higher precision than the two others, 
the average value of Op^inferred) equal to 9. 647 barns is used for all further calcula- 


TABLE HI. - VALUES OF a . OBTAINED FROM 

pot 

<7 to t CONSERVATION FROM DIFFERENT 


SOURCES OF MEASURED CROSS SECTION 


Source 


KFK-12CT 

LA-3527 b 

BNL-325 C 


°t Q t ( 20 k eV), 


13. 6 

14. 158 
13. 55 


a pot ( inferred )> 
b 

9. 500 
9.986 
9.456 


Average inferred potential scattering cross section 
>t = 9 - 647 


a Ref. 42. 
b Ref. 43. 
c Ref. 2. 
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tions (refs. 37, 38, and 39). The p-wave r n 's have been recomputed using this value 
of CT p 0 t and equation (15) to obtain the appropriate v^'s. The resulting cross section 
averages from GAROL are shown in table II under the heading DAISY (FINAL). 


Graphical Comparisons with Theoretical Distributions 

Figures 11 to 16 allow one to compare the individual resonance parameters chosen 
by DAISY to the theoretical distribution functions. These figures should indicate the 
detailed resonance structure generated by DAISY in a manner more easily grasped than 
a table of resonance parameters. The cumulative probability distributions of r n /r n 
for the s- and p-wave resonance parameters generated by DAISY are shown in figures 
11 and 12 along with the theoretical cumulative probability distributions. The cumula- 
tive probability distributions of D/D for the s- and p-wave resonance energies gener- 
ated by DAISY are shown in figures 13 and 14 along with cumulative Wigner distribution. 
It can be seen that the sampled values are in good agreement with the theoretical distri- 
butions over their range. 

While they are not comparisons of theoretical to sampled data, plots of N(E), the 
number of levels at or below the energy E, against E are shown for the s- and p-wave 
chains in figures 15 and 16. These plots serve to indicate the amount of local variation 
in level spacings. 


ANALYSIS OF DOPPLER EFFECT EXPERIMENTS: CALCULATION OF HOT TO 
COLD CAPTURE RATIOS OF THE U 238 (n,y) U 239 REACTION AT 20 keV 

The DAISY (FINAL) chain illustrated previously is used solely throughout the re- 
mainder of this report. These parameters generated by DAISY allow one to compute a 
representative shape of the cross section against energy rather than averages alone. 

In this section the parameters generated by DAISY and the code GAROL (ref. 13) will 
be used to evaluate the Doppler effect in at 20 keV. By the use of the DAISY 
(FINAL) parameters and GAROL it is possible not only to investigate the effects of reso- 
nance series overlap but it is also possible to determine separately the Doppler effect 
of the s-wave series and p-wave resonance series. 


Thin Foil Problem 

As a part of a study of the energy and temperature dependence of neutron capture 
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Figure 15. - Number of levels at or below E as a function of energy of levels for s-wave DAISY (FINAL) chain. 


measurements in the unresolved region, Seufert and Stegemann (ref. 14) measured the 

238 

hot and cold neutron capture reaction rates of a U foil, 0. 0155 centimeter thick, in 
a lead block slowing-down-time neutron energy spectrometer. The temperatures chosen 
in this experiment were 300 K for the cold sample and 750 K for the hot sample. A 
theoretical ratio of hot to cold captures is also plotted by Seufert and Stegemann as a 
function of energy from 300 eV to 30 keV. Their theoretical value of the ratio at 20 keV 
is 1. 0023. No resonance overlap corrections have been applied in their computation of 
the theoretical hot to cold capture ratio. Unfortunately, the measured ratios in this 
energy region are subject to rather large error margins. The measured ratio at about 
16. 5 keV is 1. 030±0. 024, and at about 31. 5 keV it is 1. 019±0. 045. Although experi- 
ments are very difficult to perform and experimental results are rather imprecise, it 
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is possible to compute the hot to cold reaction ratio, with some precision, in this unre- 
solved resonance region. 

The calculation of the hot (750 K) to cold (300 K) ratio for a 0. 0155-centimeter- 
thick slab of U 238 was made over the interval 19 600 to 20 400 eV by 0. 5-eV increments. 

090 

It is assumed that the slab absorber, the U foil, is embedded in a purely scattering 
medium. The results of these calculations are shown in table IV where the hot to cold 
capture ratios both with and without the overlap effect are displayed. 

It can be seen that the DAISY - GARO L calculation of the hot to cold capture ratio 
without overlap is in substantial agreement with the Seufert and Stegemann calculation 
which assumes no overlap. The value of the hot to cold capture ratio is reduced from 
1. 0030 to 1. 0023 when resonance series overlap is included. Due to the rather large 
uncertainty in the experimental results, a direct comparison between theoretical results 


27 



TABLE IV. - HOT TO COLD NEUTRON CAPTURE RATIOS 
(750 K/300 K) FOR A U 238 FOIL AT 20-keV 
NEUTRON ENERGIES 



Total cap- 
ture ratio 

R s 

R P 

With overlap 

1.0023 

±0.0016 

0.0014 

0.0009 

Without overlap 

1. 0030 
±0.0010 

.0014 

.0016 

Seufert and Stegemann theory 
without overlap 21 



1.0028 

i 


i 

Measured at 16. 5 keV a 

1.030 

±0.024 



Measured at 31. 5 keV 3- 



1.019 

±0.045 





.. J 


a Ref, 14. 


arid experiments is precluded. Yet one can determine from the detailed GAROL results 
the Doppler effects for the separate series. If one assumes that the cold capture rate 
for the s-wave resonance series is n , then the hot capture rate for the s-wave series 

o 

can be represented by n + 5n If the cold and hot capture rates for the p-wave series 

'So 

are represented in the same fashion as n and n + Sn , one can form the hot to cold 

Jr ir h* 

capture ratio as 


n s + n p + 6n s + 5n p_ 


6n„ 


= 1 + 


n s + n p 


n s + n p 


6n 


n s + n p 


5n„ 


n s + n p 


•= R, 


<5n 


P _ 


n s + n p 


R. 


( 31 ) 


( 32 ) 


The R and R values are shown for the two cases computed. It should be noted that, 
s P 

while the p-wave average capture cross section is about three times the size of the s- 
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wave average capture cross section, it (the p-wave resonance series captures) accounts 
for only about 40 percent of the Doppler effect. Assuming no overlap overestimates the 
Doppler effect for the p-wave series by about 80 percent. The reason for the rather 
small p-wave Doppler effect is that the Doppler width for a t 300 K and 20 000 eV 

is ~3. 0 eV. This is over 40 percent of the average p-wave level spacing, 7. 0 eV. 

Thus, the p-wave cross sections are much more nearly smooth than the s-wave cross 
sections. 

The error estimates shown are formed by breaking up the region from 19. 6 to 
20. 4 keV into two regions, one from 19. 6 to 20. 0 keV and the other from 20. 0 to 20. 4 
keV, computing a capture ratio for each region and then the standard deviation of these 
two capture ratios is shown below the total capture ratio. 


Thick Sphere Problem 

Another experiment utilizing a larger temperature variation and thick spheres, 

2 centimeters in diameter, has been conducted by Fieldhouse et al. (ref. 12). 

The experimental arrangement used the Li (p, n) reaction as a source of neutrons. 
The angle between the proton beam, the target, and the sample to be irradiated as well 
as the proton beam energy were chosen to obtain the required range of neutron energies 

090 

on the sample. Two samples of U were irradiated simultaneously, one sample at 
room temperature the other at 770 K. At the end of the irradiation, the hot sample was 


TABLE V. - HOT TO COLD NEUTRON CAPTURE 
RATIOS (770 K/290 K) FOR A U 238 SPHERE 
AT 20-keV NEUTRON ENERGIES 



Total cap- 
ture ratio 

R s 

R P 

With overlap 

1.023 

±0.001 

0.013 

0.010 

Without overlap 

1.030 

±0.001 

.013 

.017 

Theory of Fieldhouse et al. 
with overlap a 

1. 027 


— 

Fieldhouse et al. measured 
between 10 and 33 keV 
plotted at 21. 5 keV a 

1.034 
±0. 023 

! 

— 


a Ref. 12. 
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cooled and counted for j3 activity along with the room temperature sample. Experi- 
mental errors were reduced by alternating the temperature in the same furnaces - the 
770 K furnace is subsequently run at room temperature (290 K) and the room temper- 
ature furnace is subsequently run at 770 K. The results of these experiments are 
shown in table V. The measured value of the average Doppler effect in the energy re- 
gion, 10 to 33 keV, plotted at 21. 5 keV was 1. 034±0. 023 (ref. 12). 

The Fieldhouse theoretical results used a statistical sampling of resonance param- 
eters and a slowing down code. Therefore, this theoretical result should include s- 
wave, p-wave series overlap effects explicitly. However, rather different strength 
functions were used and a potential cross section of 10. 636 barns was used in their cal- 
culations (refs. 11 and 12). Hence, the DAISY-GAROL results could differ from the 
theoretical result of Fieldhouse et al. , a Doppler effect of 1. 027. 

The DAISY-GAROL result for the region, 19. 6 to 20. 4 keV, is 1. 023±0. 001 includ- 
ing overlap and 1. 030±0. 001 assuming no s-wave, p-wave series overlap. Thus, there 
is reasonable agreement between the DAISY-GAROL results and the theoretical result 
of Fieldhouse et al. Again the p-wave series contribution to the Doppler effect is less 
than the s-wave, p-wave cross section ratio. There is also an overestimation of the 
p-wave Doppler effect when resonance series overlap is neglected. 

The DAISY-GAROL calculations are done along the same lines as the calculations 
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of reference 12; that is, the U spheres are assumed to be covered by a thin scatter- 
ing source of neutrons. For these computations, the source is assumed to be a 0. 1- 
centimeter-thick shell around the 2-centimeter-diameter spheres. The results are not 
strongly dependent upon the thickness of this shell. When these calculations are re- 
peated for a shell with an optical thickness of 100 000. 0 centimeters, the Doppler effect, 
including overlap, is raised to only 1.046. 


SUMMARY OF RESULTS 

Using only one set of resonance parameters, determined by the code DAISY, one 
can not only match total cross sections in the unresolved energy region but one can also 
compute Doppler effects, compatible with other theoretical calculations, over a wide 
range of sample thicknesses (0. 0155 to 2 cm). 

The results obtained are as follows: 

1. The techniques used in the code DAISY allow one to efficiently choose a pseudo 
resonance parameter chain consistent with the Porter-Thomas and Wigner distribution 
functions. 

2. The cross section averages computed by the DAISY-GAROL combination are in 
very good agreement with averages computed by analytic techniques. 


30 



I 


3. A value of about 9. 6±0. 3 barns is necessary for the potential cross section in 

these calculations to preserve the measured value of the total cross section of at 

20 keV. 

4. By the use of one set of parameters, DAISY (FINAL), and the code GAROL, one 
may compute the Doppler effect at 20 keV. 

5. The Doppler effects computed by the DAISY-GAROL combination are compatible 
with other theoretical computations of the Doppler effect over a wide range of thick- 
nesses (0. 0155 to 2 cm). 

6. The effect of neglecting resonance series overlap is to overestimate the Doppler 
effect. 

7. The s-wave contribution to the Doppler effect is from 50 to 60 percent, but the 
average s-wave capture cross section is about one-third that of the p-wave capture 
cross section. 

In conclusion, a technique for the rapid generation of resonance parameter chains 

in the unresolved resonance region has been described and used in conjunction with a 

238 

standard nuclear code to evaluate Doppler effects in U 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 1, 1970, 

129-02. 
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APPENDIX A 


SAMPLING FROM A X 2 DISTRIBUTION 

The purpose of this appendix is to prove that sampling p times from a normal dis- 

2 

tribution, squaring the samples and summing provides a x sample of p degrees of 
freedom. 

First, M x (0), the moment generating function of f(x) is defined as follows: 


Mjd) = 



e^ x f(x)dx 


Also, it can be shown (refs. 29, 30, and 40) that, if two random variables have the same 
moment generating function, they are identical. 

o 

Next, the moment generating function of the x distribution with p degrees of 
freedom shall be computed as follows: 


M 2 (d) 
X 




e 0X 2 (x 2 ) ( I ,-2)/2 e -x 2 /2 


d(x 2 ) 


If 


is replaced by w, the result is 


M w («) = 



/" e 0w>-2)/2 e -3/2 dw 

0 


If z = w(l - 26>)/2 dw = 2 dz/(l - 20), 


M wW = 11^ /"e-W 2 )-^ 

F © ° 


M^(e) = t 1 - - 2| » " /2 t(a\ 

r (i) W 
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M w (e) = (i - 2 ey v / 2 
M n(0) = (i - 2 e) - "/ 2 

x 2 

Now the moment generating function of the sum of the squares of samples chosen 
from the normal distribution function will be determined. Use will be made of the fol- 
lowing theorem (refs. 29, 30, and 40). The moment generating function of the sum of 
n independent random variables is equal to the product of the moment generating func- 
tions of the individual variables. If x is normally distributed with mean 0 and vari- 
ance 1, let a random sample of v values of x^ be taken from the distribution. Now 
the moment generating function of q shall be computed. Where 


q = Z xf 

i=l 


M (0) = M 2 2 ( 0 ) = M 2 (0) - M 2 (0) . . . M 2 (0) = [m 2 (e)l" 

4 Xj + . . . +X X 1 Xo X L X J 


Since x is normally distributed, 

m 2 {e)=^— f 
yr y/2n J_ x 

Let y = x(l - 2 O)^ 2 ; then 


e** 2 e- 2 / 2 dx 


M 2 (0) = (1 - 20)” 1 / 2 J— f e‘ y2//2 dy = (1 - 20)' 1 / 2 
yr V2 * 


Recalling the previous result gives 


M q (0) = = [(1 - 26)~ 1/2 J = (1 - 20)' 


v/2 


2 

which is exactly the same result as the generating function of x. with v degrees of 
freedom. 
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Thus, since the moment generating functions are equal, the distribution functions 

2 

are identical. If one wishes to sample from a x distribution with v degrees of free- 
dom, one needs only to sample the normal distribution v times, square the results, 
and sum. 
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APPENDIX B 


FAST GENERATOR FOR NORMAL RANDOM VARIABLES 


2 

It has been shown that, if one wishes to sample from the x distribution with p 
degrees of freedom, the problem may be reduced to sampling from the normal distribu- 
tion. The purpose of this appendix is to outline in some detail a fast procedure (ref. 23) 
for sampling the normal distribution. 

The procedure is most simply explained as follows. The density function f(x) of 
the positive normal distribution is divided into three portions 


f(x) = 


V2 e- 1 / 2 * 
Vtt 


= 0. 9578 f X (x) + 0. 0395 f 2 (x) + 0. 0027 fg(x) 


The meaning of this formula is that one generates, with probability 0. 9578 a random 
variable with density f^, with probability 0. 0395 a random variable with density f 2 , 
and finally with probability 0. 0027 a random variable with density fg. The exact prob- 
abilities can be determined from reference 23. The importance of this splitting of the 
density function is that the procedure which gives f ^ is extremely fast, the procedure 
for f 2 is rather short, and the time for fg rather long. 

The procedure for f^ depends upon representing much of the area under the den- 
sity curve f(x) as a series of uniformly distributed random variables with differing 
probabilities of being chosen (sketch (a)). A stored series of values then allows the use 
of the first few digits of a uniformly ‘distributed random number to choose the correct 
leading digits of the variate f^ and the rest of the uniform random variable gives the 
rest of f^. If the stored values are contained in the vector A i and if the uniformly 
distributed random variable u has digits 0 • u^, u 2 , Ug, . . . , u^. If 00 < u i u 2 — ?9, 

put the normal random variable x = A„ „ + ■ 0 __ „ 

u l u 2 U 4 U 5 U 6 U 7* 

If 790 < u-UoUq < 940, put x = A - 771 + • 0„ „ „ „ . This will generate the 

1 z 6 U 1 U 2 U 3 U 4 U 5 6 7 

variable shown in sketch (a). In order to generate the correct random variable f(x) in 

the range 0 < x < 3 . 0 , a sample must also be taken from the density function f 2 (x) 


1.0 

.8 

f l <x) *4 
0 


TTriT n-r^ . i 

1.0 1.5 2.0 2.5 3.0 


(a) 
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0 .5 1.0 1.5 2.0 2.5 3.0 


x 

(b) 






(c) 


shown in sketch (b) with probability 0. 0395. This is accomplished by a modified rejec- 
tion technique. It can be shown for a nearly linear density function, as is shown in 
sketch (c), that the following technique will provide a correct sampling (ref. 41). In 
order to generate a random variable Y from a nearly linear density function g(x) for 
t < x < t + a, first contain g(x) within two parallel lines. Then choose two independent 
random variables u and v distributed uniformly on the unit interval (0, 1). If the 
larger of u and v is less than If. h, set Y = t + a min(u, v). If the larger of u and v 
is not less than l/h, then find out if 

h |u - v | < g[t + a min(u, v)] 

If true, put Y = t + a min(u, v); otherwise generate a new u and v and start this en- 
tire procedure again. 

The two techniques described previously allow one to generate the correct density 
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function for the normal distribution to a value of the random variable of 3. 0. Beyond 
3.0 a process which is executed for only 0. 27 percent of the values. It is a modification 
of a technique found in reference 25 and shall not be described inasmuch as it does not 
affect the timing a great deal. The interested reader is referred to references 23 
and 25. 

The aforementioned procedures including the one not described in detail have been 
programmed in a FORTRAN IV subroutine. It provides on the order of 5000 normal 
random variables per second on the Lewis 7094-11. 
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